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In this paper we give a method to associate a graph with an arbitrary density ma- 
trix referred to a standard orthonormal basis in the Hilbert space of a finite dimensional 
quantum system. We study the related issues like classification of pure and mixed states, 
Von-Neumann entropy, separability of multipartite quantum states and quantum oper- 
. . , ations in terms of the graphs associated with quantum states. In order to address the 

^ ! separability and entanglement questions using graphs, we introduce a modified tensor 

O I product of weighted graphs, and establish its algebraic properties. In particular, we show 

that Werner's definition [1] of a separable state can be written in terms of a graphs, for 
^ I the states in a real or complex Hilbert space. We generalize the separability criterion 

(degree criterion) due to S.L. Braunstein, S. Ghosh, T. Mansour, S. Severini, R.C. Wilson 
[2], to a class of weighted graphs with real weights. We have given some criteria for the 
^ I Laplacian associated with a weighted graph to be positive semidefinite. 

> : PACS numbers:03.67.-a,03.67.Mn 

m . 
in ; 

(N ' 1 Introduction 

o . 
\o . 

■ Quantum information is a rapidly expanding field of research because of its theoreti- 
^ ■ cal advances in fast algorithms, superdence quantum coding, quantum error correction, 

teleportation, cryptography and so forth [3,4,5]. Most of these applications are based 
. on entanglement in quantum states. Although entanglement in pure state systems is 

^ ■ relatively well understood, its understanding in the so called mixed quantum states [6], 

O^' which are statistical mixtures of pure quantum states, is at a primitive level. Recently, 

■ Braunstein, Ghosh and Severini [2, 7], have initiated a new approach towards the mixed 
^ ■ state entanglement by associating graphs with density matrices and understanding their 

■ classification using these graphs. Hildebrand, Mancini and Severini [8] testified that the 
degree condition is equivalent to the PPT-criterion. They also considered the concurrence 
of density matrices of graphs and pointed out that there are examples on four vertices 
whose concurrence is a rational number. In this paper we generalize the work of these au- 
thors and give a method to associate a graph with the density matrix (real or complex), of 
an arbitrary density operator, and also to associate a graph with the matrix representing 
hermitian operator (observable) of the quantum system, both with respect to a standard 
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orthonormal basis in Hilbcrt space. We define a modified tensor product of graphs and 
use it to give Werner's definition for the separabihty of a m-partite quantum system, in 
M'?i (g)M''2 (g, . . . (gjRgm^ as well as ^C^^ ^ . . .(^£<im terms of graphs. We also deal with 
classification of pure and mixed states and related concepts like Von- Neumann entropy in 
terms of graphs. 

The paper is organized as follows. In Section 2, we define weighted graphs and their 
generalized Laplacians which correspond to density matrices, and discuss the permutation 
invariance of this association. We also deal with pure and mixed states in terms of graphs. 
Section 3 deals with Von- Neumann entropy. Section 4 is concerned with separability issues 
as mentioned above. In Section 5, we deal with graph operations which correspond to 
quantum operations [5, 9, 10]. In Section 6, we present a method to associate a graph 
with a general hermitian matrix, having complex off-diagonal elements. We define the 
modified tensor product for complex weighted graphs and express the separability of mixed 
quantum states in a complex Hilbert space in terms of graphs, using Werner's definition. 
In section 7, we present some graphical criteria for the associated Laplacian to be positive 
semidefinite. Finally, we close with a summary and some general comments. Sections 2 
to 5 deal with graphs with real weights , that is, quantum states living in real Hilbert 
space. Graphs with complex weights, corresponding to density operators with complex off 
diagonal elements are treated in section 6. However, a large part of the results obtained 
for real Hilbert space in sections 2 to 5, go over to the case of complex Hilbert space (see 
section 8 (ix)). 

2 Density matrix of a weighted graph 
2.1 Definitions 

A graph G = {V, E) is a pair of a nonempty and finite set called vertex set V{G), whose 
elements are called vertices and a set E{G) C V'^{G) of unordered pairs of vertices called 
edges. A loop is an edge of the form {vi, Vi} for some vertex Vi. A graph G is on n vertices 
if \V{G) I = n. We call the graph as defined above a simple graph. \E{G) \ = m + s, where 
m : number of edges joining vertices, s : number of loops in G [11]. 

A weighted graph {G, a) is a graph together with a weight function [12] 

a : V{G) X V{G) R 

which associates a real number (weight) a{{u, v}) to each pair {u, v} of vertices. The 
function a satisfies the following properties: 

(i) a{{u, v}) if {u, v} e E{G, a) and a{{u, v}) = if {u, v} ^ E{G, a). 



2 



(ii) a{{u,v}) = a{{v,u}) 



(iii) a{v, v) if {v, v} e E{G, a) and is zero otherwise. 

If e = {u,v} is an edge in E{G,a), property (ii) allows us to write a(e) or auv for 
a{{u,v}). A simple graph can be viewed as a weighted graph with all nonzero weights 
equal to 1. 

In the case of simple graphs the degree of a vertex v e V{G) is defined as the 
number of edges in E{G) incident on v. For a weighted graph we set 

d{G,a){v) ^ ^ ^ a„„. (1) 

ueV{G,a) 

The adjacency matrix of a weighted graph with n vertices M{G,a) = [a,uv]u,vev(G,a) is 
a, n X n matrix whose rows and columns are indexed by vertices in V{G, a) and whose 
uv-th element is a^v Obviously the adjacency matrix M{G, a) is a real symmetric matrix 
with diagonal element vv equal to the weight of the loops on vertex v (i.e a^^). 

The degree matrix for the weighted graph A(G, a) is a n x n diagonal matrix, whose 
rows and columns are labeled by vertices in V{G, a) and whose diagonal elements are the 
degrees of the corresponding vertices. 

A{G,a) ^diag[d^;v eV{G,a)]. (2) 

The combinatorial Laplacian of a weighted graph is defined to be 

L{G,a)^A{G,a)-M{G,a). (3) 

The degree sum of (G, a) is defined as 

diG,a)= d,^TrA{G,a). (4) 

veV{G,a) 

The Laplacian defined by (3) has no record of loops in the graph. Therefore we define 
the generalized Laplacian of a graph {G, a), which includes loops, as 

Q{G,a) = A{G,a) - M{G,a) + Ao{G,a) (5) 

where Ao(G, a) is a n x n diagonal matrix with diagonal elements equal to the weights of 
the loops on the corresponding vertices 

[Ao(G,a)]„„ = a„„. (6) 

We call Ao(G, a) the loop matrix of the graph {G, a) . 
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For a given weighted graph {G,a), the generahzed Laplacian, defined by (5), is not 
necessarily a positive semidefinite matrix. When, for a given graph {G,a), the general- 
ized Laplacian Q{G, a) is positive semidefinite, we can define the density matrix of the 
corresponding graph (G, a) as 

a{G, a) = -^Q{G, a) = ^[L{G, a) + Ao{G, a)] (7) 

where Tr{a{G, a)) — 1. Note that, this definition of the density matrix of a weighted 
graph (G, a) reduces to that of the density matrix for a simple graph without loops [7] . 

Whenever we can define the density matrix for a graph (G, a) we say that the graph 
(G, a) has density matrix. 

For any density matrix cr, we can obtain the corresponding graph as follows: 

(i) Determine the number of vertices of the graph from the size {n x n) of the density 

matrix. The number of vertices = n. Label the vertices from 1 to n. 

(ii) If the ij-th element of cr is not zero draw an edge between vertices Vi and Vj with 
weight — CTjj. 

(iii) Ensure that dy^ — an by adding loop of appropriate weight to Vi if necessary. 
Example{l): For the following three matrices, we find the corresponding graphs. 

(l)a = 



1 1 
1 1 



in 
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Figure 2 
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Figure 3 

2.2 Invar iance under isomorphism 

Two weighted graphs (G, a) and (G", a') are isomorphic if there is a bijective map [13] 

: V{G,a) I — > V{G',a') 

such that 

{cl){vi), (l>{vj)} e £;(G"a') iff {Vi, vj} e a), i, i = 1, 2, • • • , n 

and 

We denote isomorphism of {G.a) and (G'.a') by (G, a) = {G'a'). 

Equivalcntly, two graphs (G, a) and (G', a') are isomorphic if there exists a permutation 
matrix P such that 

P^M(G, a)P ^ M(G',a'). 

Note that, 

P^A(G, a)P = A(G', a'); P^Ao(G, a)P = Ao(G', a') 
Therefore we have 

P^Q{G,a)P = Q{G',a'). (8) 
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This means that Q{G, a) and Q{G', a') arc similar and have the same eigenvalues. There- 
fore, if Q{G,a) is positive semidefinite then so is Q{G',a'). Therefore, if {G,a) has the 
density matrix so does (G", a'). We have proved 

Theorem 2.1 : The set of all weighted graphs having density matrix is closed under 
isomorphism. ■ 
Since isomorphism is an equivalence relation, this set is partitioned by it, mutually 
isomorphic graphs forming the partition. 

2.3 Correspondence with quantum mechanics 

Henceforth, we consider only the graphs having density matrix unless stated otherwise. 
The basic correspondence with QM is defined by the density matrix of the graph. For a 
graph with n vertices the dimension of the Hilbert space of the corresponding quantum 
system is n. To establish the required correspondence we fix an orthonormal basis in 
the Hilbert space C?) CS> • • • M'*" of the system, which we call the standard basis 
and denote it by {\ijM ■ ■ j, k, i ■ ■ ■ = 1,2, ■■■ ,n = qiq2---qm, or by {\vi)},i = 
1, ■ ■ ■ ,n = qiq2 ■ ■ ■ qm- We label n vertices of the graph (G, a) corresponding to the given 
density matrix by the n basis vectors in the standard basis. We say that the graph {G, a) 
corresponds to the quantum state (density operator) whose matrix in the standard basis 
is the given density matrix. Finally, we set up a procedure, by associating appropriate 
projection operators with edges and loops of {G, a) to reconstruct this quantum state 
from the graph (G, a). (See Theorem 2.7). In view of Theorem 2.1 , if (G, a) has density 
matrix a and {G, a) = (G", a') with the corresponding permutation matrix P,then {G', a') 
has the density matrix P^aP. All of this paragraph applies to the complex weighted 
graph (section 6). 

2.3.1 Pure and mixed states 

A density matrix p is said to be pure if Tr{p'^) = 1 and mixed otherwise. Theorem 
2.3 gives a necessary and sufficient condition on a graph {G, a) for cr{G, a) to be pure. 
For a graph (G, a) having k components (Gi, ai), (G2, 02)1 ■ ■ ■ i (Gk, dk) we write (G, a) = 
(Gi, ai) l±) (G2, 02) tbi ■ ■ ■ l±l (Gfc, ttk) where ai,i = 1, ■ ■ ■ , k are the restrictions of the weight 
function of the graph (G, a) to the components. We can order the vertices such that 
M(G, a) — ®|^iM(Gi, Oj). When k — 1, (G, a) is said to be connected. From now on we 
denote by Xi{A), A2(^), • • • , Afc(A) the k different eigenvalues of the Hermitian matrix A in 
the nondecreasing order. The set of the eigenvalues of A together with their multiphcities 
is called the spectrum of A [13, 14, 15]. 

Lemma 2.2 : The density matrix of a graph (G, a) without loops has zero eigenvalue with 
multiplicity greater than or equal to the number of components of (G, a) with equality 
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applying when weight function a =constant > 0, 

Proof : Let (G, a) be a graph with n vertices and m edges. Since Q{G,a) is positive 
semidefinite, for x e iR" we must have [12] 

m s 

x^Q{G, a)x = ai^j^{xi^ - Xj^f + ^ 0^,^,4 ^ 0- 
k=i t=i 

For the graph without loops the above inequality becomes 

x^QiG, a)x = tti^j^ixi^ - Xj^y > 0. (9) 

For x'^ = {1 1 • • • 1) we can see x^Qx — 0. This means that x'^ — {1 1 1 • • • 1) is an 

unnormalizcd eigenvector belonging to the eigenvalue [13]. If there are two components 
(Gi,a) and {G2-(i) of (G^a)^ with ni,mi and 722,7712 as the number of vertices and edges 
in {Gi,a) and {G2,a) respectively, we can decompose the sum in equation (9) as 

mi m2 

x'^Q{G,a)x= XI "^/^lifei (^ifei - ^j^f + ^ (^ik^jk,i^ik, -^jk,f- (10) 

fcl=l fc2=l 

For a:^ = (1 1 1 • • • 1) both the terms in (10) vanish. Now consider two vectors xj — 
(0 • • • 01 1 • • • 1) with first rii components zero and last n2 components 1 and = 
(1 1 • • • 1 • • • 0) with first rii components 1 and last n2 components zero, (721+712 = n). 
Obviously the RHS of (10) vanishes for both xi and X2. This imphes xi and X2 are two 
orthogonal eigenvectors with eigenvalue zero. This means multiplicity of zero eigen value 
is at least 2 (number of components in (G, a)). 

The equality condition for a^^ =constant> V {u,v} G E{G.a) is proved in [7]. ■ 
Theorem 2.3 : The necessary and sufficient condition for the state given by a graph 
{G, a) to be pure is 

n m 

Y.dl + ^T.<j.-dka) (11) 

1=1 fc=l 

where di is the degree of the vertex Vi , ai^j,^ is the weight of the edge {vi^,Vj^}, (vj^. 7^ vj^) 

n 

and (i(G,a) is the degree sum d(G,a) = ^2^1- 

i=l 

Proof : Equation (11) is just the restatement of the requirement Tr{a'^{G, a)) — 1, which 
is the necessary and sufficient condition for the state a{G, a) to be pure. ■ 
Lemma 2.4 : The graph {G, a) for a pure state a{G, a) has the form {K^, b) l±) v^+i l±) 
^^+2 l±) • • • l±) for some 1 < I <n. 
Proof : Since the state is pure, it has the form 

(. 

1^) = ^(^ik\'^ik)^^ <ik<n. 

k=l 
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We can permute the basis vectors to transform this sum to = '^Ci\vi). That is, 

i=l 

the i basis kets contributing to the sum in the above equation become the vectors 
ki))k2),-'" ) l'^^^) under this permutation. The resulting density matrix has a 

block of first i rows and first i columns all of whose elements are nonzero, while all the 
other elements of density matrix are zero. The graph corresponding to this density matrix 
is just the required graph. ■ 
Example (2) : We now give important cases of pure state graphs in which we use 
later. 



Figure 4 
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sponding graph is as shown in Figure 4. 
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P[|t'i)], the corresponding graph is as shown 
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The corresponding graph is as shown in Figure 6. 



^'[75(K) + I^2))],a>0. 




Figure 7 
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(iv) 



<7(G,a) = - 



1-1-11 

-11 1-1 

-11 1-1 
1-1-11 



-P[{\-)\-))] 

, where |— ) = ~ 1 2)), the corresponding graph is as shown in Figure 7. 




Figure 8 
(v) 



a{G,a) 



1 
4 



110 110 

110 110 



110 110 

110 110 





where |+) = + in ® 1^^, the corresponding graph is as shown in Figure 8. 

It may be seen that in each of the cases in example (2), same density matrix on the 
standard basis corresponds to infinite family of graphs as the nonzero weight on each edge 
or loop is multiplied by a constant. But this is a false alarm because any weight a 7^ 1 
only changes the length of the corresponding state vector in the Hilbert space (i.e. state 
becomes unnormalized) which does not have any physical significance. Another example 
pertaining to this situation is the random mixture (see Lemma (3.1 )). 



a(G, a) = — 
an 
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-In 

n 



However, this does not lead to any contradiction because of the uniqueness of the random 
mixture [6]. 
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All the density matrices in (i), (ii), (iii), (iv), (v) above represent pure states. 
Remark 2.5 : Any graph with weight function a = constant > has the same density 
matrix for all a > 0. This infinite family of graphs corresponds to the same quantum 
state (density operator). 

Definition 2.6 : A graph {H, h) is said to be a factor of graph [G, a) if V{H, h) = V{G, a) 
and there exists a graph (//', h') such that V{H', h') = V{G, a) and M{G, a) = M{H, h) + 
M{H',b'). Thus a factor is only a spanning subgraph. Note that 



hviVj if {vhVj] e E{H,b) 
b' if {v.,v,}eE{H',b') 



Now let (G, a) be a graph on n vertices Vi, ■ ■ ■ , t'„ having m edges 
{vi^,Vj^}, ■■■ , {vi^,Vj^} and s loops {vi^,Vi^} ■ ■ ■ {vi^,Vi^} where 1 < iiji, • • • , imjm < 
n,l < iii-2 ■ ■ - is ^ n. 

Let {Hi^ji^, aii^j^) be the factor of {G, a) such that 

\M(H n )^ -/ ^'^^^ ""^'^ t^^jk or u^jk,w^ik . . 

Let {Hi^^i^, tti^i^) be a factor of {G, a) such that 

\M(H n ^] _\ (^kk when u^it^w , . 

[M{H,,,,,a,,,,)U-^ otherwise 

Theorem 2.7 : The density matrix of a graph (G, a) as defined above with factors given 
by equation (12) and (13) can be decomposed as 

^ m 1 * 

(^{G, a) = 2ai^j^a{Hi^j^, Ui^jJ + ai^i^a{Hi^i^, a^^jj (14) 

«(G,a) ^ a(G,a) 

or 

^ m ^ 1 * 

a{G,a) = - J]2a,,,,P[-=(|^,J - \vj^))] + - J] a,,,,P[|^,,)] (15) 

Proof : From equation (12), (13) and Theorem 2.3 and Lemma 2.4, the density matrix 

is a pure state. Also, 

(^{Hitk-,0'kk) — [^o{Hit,k->0^kk)] 

O'kk 



10 



is a pure state. Now 



k=l t=l 
m s 

M(G, a)^Yl ^(HikJk^ + Yl MHuu,ai,i,). 



k=l 



t=l 



Therefore 



a{G, a) 



(G,a) 



.k=l 



k=l 



+ 



(G,a) 



t=\ 



^ m 1 * 

d{G,a) ^ Ct(G,a) 



1 



(14) 



J / ^ ^kJk" \''i-kJkT "i-kUkJ • J 

d(G,a) ^ d^G,a) ^ 

In terms of the standard basis, the uw-th element of matrices a{H.i^j^,ai^j^) and 
a{Hi^i^,ai^i,) are given by {vu\cr{Hi^j^, , ai^j^)\vyj) and {vu\cr{Hi^i^ai^i^\vyj) respectively. In 
this basis ^ 

'^iHiuk,aikjk) = P[^iM - \^jk))] 

Therefore equation (14) becomes 

1 ™ 1 1 * 
a{G, a) = 2«iu.^[^(l^iJ - K)) + E ^hhP\W)\ (15) 



Remark 2.8 : If all weights di^j^, > then equations (14), (15) give a{G, a) as a mixture 
of pure states. However, in the next subsection we show that any graph {G, a) having 
density matrix can be decomposed into graphs (spanning subgraphs) corresponding to 
pure states. 



2.3.2 Convex combination of density matrices 

Consider two graphs {Gi,ai) and (G'2-Q'2) each on the same n vertices, having a{Gi,ai) 
and <7{G2,a2) as their density matrices respectively. We give an algorithm to construct 
the graph (G, a) whose density matrix is 

a{G, a) = Xa{Gi, oi) + (1 - X)a{G2, as) 
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0<X<1,\ = a /l3, a, f3>0 are real. 

We use the symbol U to denote the union of the edge sets of two graphs (Gi, ai) and 
(^2,02) on the same set of vertices to give (G.a). If {vi,Vj} G E{Gi,ai) and {vi,Vj} G 
E{G2,a2) then a{{vi,Vj}) = ai{{vi,Vj}) + a2{{vi,Vj}). We write (G, a) = (Gi,ai) U 
(^2,02)- If E{Gi,ai) and E{G2,a2) are disjoint sets, then we call the resulting graph 
{G, a) the disjoint edge union of (Gi, ai) and {G2, 02), we write {G, a) — {Gi, ai)-i-(G2, 02). 

The algorithm is as follows : 
Algorithm 2.9 : 

(i) Put X = a/(3 so that (1 — A) = where a > ,(3 > are real. 

(ii) Write a{G, a) = ^{aa{Gi, oi) + {P - a)a{G2, 02)). 

(iii) Modify the weight functions of the two graphs (Gi, ai) and (G2, 02) to get a[ = aai 
and 02 = (/? — a)a2. 



(iv) The graph (G, a) corresponding to a in step (ii) is 

(G,a) = (Gi,a;) U (Gs,^) 

such that 



(16) 



(Ol)^.^,. + {(l2)viVj 



(16a) 

a^,^, = (a'J^i^^i + W2)viVi (166) 
where we take {a[^2)viVj = = (a'i_2)<^ii'i if {vi,Vj}, {vi,Vi} ^ E{Gi,ai) or E{G2,a2) 



We can apply this algorithm to any convex combination of more than two density 

k 

matrices a{G, a) — '^PiCr{Gi, ai), — 1, by writing pi — ai/(3, ai, (3 > and real, i — 

1=1 i 

l,---,k. 

Example(3) : consider the density matrices 



1111 
1111 
1111 
1111 
whose graph is shown in Figure 9 



(i) (j(Gi,ai) = ! + +)(+ + 



12 



Figure 9 



(7(^2,02) = 



whose graph is shown in Figure 10 
11 12 




21 

Figure 10 



l|ii)(ii| + i|V^+)(V^^ 



1 

4 



2 

110 

110 





The graph corresponding to 

a{G, a) 



1 2 
-a{Gi,ai) + -a{G2,a2) 



1 

12 



5 111 

13 3 1 
13 3 1 
1111 



is given in Figure 11 
Lemma 2.10 : Let (Gi,ai), ((^2,02) and (G, a) satisfy 



or, 



{G,a) = (Gi,ai) U (^2,02) 



(G,a) = (6-1,01) + (62,02). 



13 



Figure 11 



Then 

Q{G,a) = Q{Gi,ai) + Q{G2,a2) 

and 

ct(G, a) = — oi) + — (t(G2, 02). 

"{G,a) a(G,a) 

Proof : For two factors of (G, a), (Gi, oi) and (6*2, 02) we have 

M{G,a) = M{Gi,ai) + M{G2,a2) 

A{G,a) = A{Gi,ai) + A{G2,a2) 

Ao(G,a) = Ao(Gi,ai) + Ao(G2,a2) 

L{G,a) = A{G,a)-M{G,a) 

Q(G,a) = L(G,a) + Ao(G,a) 

Substitute M{G, a), A{G, a), Ao(G, a) and L{G, a) in Q{G, a) as above to get 

Q{G,a)^Q{G,,a,) + Q{G2,a2) 

and also 

o[G, a) = a(Gi, ai) + — o-(G2, 02). 

"{G,a) a(G,a) 

■ 

Remcirk 2.11 : Obviously, the operation U is associative. We can apply Lemma 2.10 for 
more than two graphs, 

(G, a) = Ui{Gi, tti) ^ Q{G, a) = ^ Q{Gi, a^) 

i 

and ^ 

a{G, a) = V! ^(^i' 

«(G,a) ^ 

Theorem 2.12 : Every graph (G, a) having density matrix a{G,a) can be decomposed 
as (G, a) = Uj(Gi, a^) where (T(Gi, a^) is a pure state. 
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Proof : Every density matrix can be written as the convex combination of pure states 

k 

a(G,a) = Y.PMi'^il- 

By applying Algorithm (2.9), Lemma 2.10 and Remark 2.11, we get the result. ■ 

2.3.3 Tracing out a part 

Consider a bipartite system with dimension pq. Let a{G, a) be a state of the system 
with graph (G, a) having pq vertices labeled by — 1, • • • ,p and j — 1, • ■ ■ ,q. li we 

trace out the second part with dimension q, we get the state of the first part which is 
p X p reduced density matrix of o-{G, a). The corresponding graph {G', a') has p vertices 
indexed by (i) and its weight function a' is given by 

9 



k=l 



and 



k=l leV{G',a') 

Where dik is the degree of vertex (ik) in original graph. 

Example (4) ■ Consider a graph (G, a) as shown in Figure 12a in 
corresponding density matrix is 



e. The 



<7^^(G,a) = - 



9-1-11 

-1 3 -1 -1 

-1 -1 3 -1 
1-1-11 



after tracing out the second particle the graph {G',a') on two vertices becomes as in 
Figure 12b with corresponding density matrix 



1 


"12 -2 " 


1 


"6 -1 " 


16 


-2 4 


~ 8 


-1 2 



a^{G',a') = - 
which is the same as the reduced density matrix of a^^. 

3 Von Neumann entropy 

The Von Neumann entropy of n x n density matrix a is 

n 



15 




a 

Figure 12a 



Figure 12b 

It is conventional to define OlogO = 0. The Von Neumann entropy is a measure of 
mixedness of the density matrix. For a pure state cr, S{a) — 0. 



3.1 Maximum and minimum 

Let 

{G,a) = li^U{Kla,) (17) 

where {Kl, a,j) is the graph on i-th vertex with a loop having weight > 0. 
Lemma 3.1 : Let (G, a) be given by (17) with the additional constraint that ai = c — 
— 1,2,-- - ,n. The density matrix of the graph (G, a) is the random mixture of pure 
states with a{G, a) — ^In- 

Proof : For the graph (G, a), the first term in (14) vanishes. Then 



1 " 

'^{G, a) = Ao(-f/"i,i,, a) 

d{G,a) ^ 



where Ao(i?j(i(, a) is the nxn matrix with all elements zero except {it, it)th element which 
is equal to a. This means 



a{G,a) 



d{G,a) 







n 



because d, 



na. 



Theorem 3.2 : Let (G, a) be a graph on n vertices. Then 
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(i) max(G,a) S{a{G, a)) = logsn 

(ii) min((3 a) s{a{G, a)) — 0, and this value is attained if a{G, a) is pure. 

Proof : (i) By Lemma 3.1 a{G,a) defined in the Lemma has eigenvalue 1/n with multi- 
plicity n. The corresponding Von Neumann entropy is log2 n. Since {G, a) is on n vertices, 
the support of a{G, a) has dimension < n. Any matrix having dimension of support < n 
cannot have Von Neumann entropy > log2 n. 

(ii) For pure state S{a{G, a)) = and S{a{G, a)) ^ . ■ 

4 Separability 

In this section we primarily deal with the graphs representing a bipartite quantum system 
with Hilbert space R^^IR'^ of dimension pq. Obviously, the corresponding graph has n = pq 
vertices. We label the vertices using standard (product) basis {\vi) = \us+i) ® \wt)},0 < 
s<p — l,l<t<q,i = sq + t. 

4.1 Tensor product of weighted graphs 

The tensor product of two graphs {G, a) and {H, h) denoted (G, a) ® {H, h) is defined as 
follows. 

The vertex set of (G, a) ® {,H,h) is V{G,a) x V{H,b). Two vertices {ui,Vi) and 
(^2, V2) are adjacent if {^i, U2} G E{G, a) and {vi, V2} G E{H, h). The weight of the edge 
{{ui,Vi),{u2,V2)} given by a{„i,^2}6{^,,„2} and is denoted by c({(iii, Vi), (it2, ^2)}). Note 
that either ui and U2 or vi and V2 or both can be identical, to include loops. 

The adjacency, degree and the loops matrices of (G, a) ® {H, h) are given by 

M((G, a) ® {H, b)) = M{G, a) ® M{H, b) (18a) 

A((G, a) (g) {H, h)) = A(G, a) ® A{H, b) {18b) 
Ao((G, a) (g) (//, b)) = Ao(G, a) Ao(//, b) (18c) 

Note that 

L((G, a) ® (H, b)) ^ L{G, a) (g) L{H, b) 

Q{{G, a) ® {H, b)) ^ Q{G, a) ® QiH, b). 

In fact, in general, the tensor product of two graphs having density matrix may not have 
density matrix. 

For two simple graphs G and H we know that [16, 7] 

da^H — da ■ dn- 
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This result is also satisfied by the tensor product of the weighted graphs. 

d{G,a)®{H,h) = d(G,a) ' d(H,b)- (19) 

4.2 Modified tensor product 

We modify the tensor product of graphs in order to preserve positivity of the generalized 
Laplacian of the resuhing graph. 

Given a graph (G, a) we define {C^, a) by 

^(6"^, a)^V{G,a) 

E[{G^, a)] = E{G, a) \ {{v,, v,} : {v,, v,} G E{G, a)} 

That is (C^, a) is obtained from {G, a) by removing all loops. 
Given a graph {G, a) we define {G, a) by 

V{G,a)^V{G,a) 

E(G, a) = E{G, a) \ {{vi, Vj} : i ^ j, {v„ Vj} e E{G, a)}. 

That is, {G, a) is obtained by removing all edges connecting neighbors and keeping loops. 

Note that in both {G'^,a) and {G,a), weight function a remains the same, only its 
support is restricted. 

Given a graph (G, a) we define (— G, a) = (G, — a). Given a graph (G, a) we define 
(G#,a') 

V{G*,a') = V{G,a) 
{G*,a') = ^-{K,,a',) 

where Ki is the graph consisting of ith vertex with a loop and a'^ is the weight of the loop 
on the ith vertex. If = then there is no loop on the ith vertex. a'^,i — 1,2, ■■ -n are 
given by 

VkeV{G,a) 

Note that the term Vk — Vi is also included in the sum. 
We now define the graph operators on the set of graphs 



(i) T] : {G,a) ^ {-G,a) = {G,-a) 

(ii) C: {G,a) ^ {G'^,a) 

(iii) AA : (G, a) (G#, a') 

(iv) n : (G, a) (G, a) 



(206) 
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a2 



^ (a) • ■ ^^(b) • 



1 , 2 ^^^^ 1 2 ^2 1 2 



(d) 



(c) 

Figure 13 



Some properties of the graph operators defined in (20b) are, 



(i) M{r]{G,a)) = -M{G,a) 

A{ri{G,a)) = -A{G,a) (21) 
Ao(r?(G,a)) = -Ao(G,a) 



(22) 



(23) 



dr,{G,a) = -d{G,a) 

(ii) M{/:{G,a))^M{G,a)-Ao{G,a) 
A{C{G,a))^A{G,a)-Ao{G,a) 

Ao(>C(G,a)) = [0] 

c?/;(G,a) = Tr{A{G, a)) - Tr{Ao{G, a)) = d^c^a) 

(iii) M{Ar{G,a)) = A{G,a) 
A{X{G,a)) = A{G,a) 
Ao{J\f{G,a))^A{G,a) 
dM(G,a) = Tr{A{G, a)) = d^G,a) 

(iv) M(Q(G',a)) = Ao(G',a) 

A((](G,a)) = Ao(G,a) 
Ao(fi(G',a)) = Ao(G,a) 
c?n{G,a) = Tr(Ao(G',a)) 

Example (5) : Given a graph {G, a) as shown in Figure 13a, if we act by rj, C ,J\f and 
Q on (G, a), wc get the graphs ri{G, a), C{G, a), NiG, a) and Vt{G, a) as shown in Figures 
13b, 13c, 13d and 13e respectively. 

Definition 4.1 : Let (G, a) and (if, h) be two graphs with p and g(> p) vertices respec- 
tively. Then their modified tensor product is defined by 

(G,a)H(i/,6) = {C{G,a)® Cr]{H,h)} + {C{G,a)®J\f{H,h)] 

+{AA(G, a) ® C{H, b)} + {fl{G, a) ^ n{H, b)} (25) 



(24) 
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V{{G, a) B {H, b)} = V{G, a) x V{H, b) 

whose cardinality is pq. 

E{(G, a) □ {H, b)} = Disjoint union of the edge set of each term in (25). 
Lemma 4.2 : (i) A{{G, a) B {H, b)) = A{G, a) ® A{H, b). 

(ii) Ao((G, a) B (H, b)) = Ao(G, a) ® Ao(i/, b). 
Proof : Consider the degree matrix of the modified tensor product we have 

A{{G,a)B{H,b)) = A{jC{G,a)^jCr){H,b)) + A{jC{G,a)^Af{H,b)) 

+A{Af{G, a) (g) jC{H, b)) + A{n{G, a) n{H, b)) 



This follows from Lemma 2.10. Using equation (18b) and equations (21) to (24) to the 
terms on the RHS of the above equation we get 

A{{G, a) B {H, b)) = A{G, a) ® A{H, b). 

Equations (ii) can be proved similarly. ■ 
Corollary 4.3 : c?(G,a)[3(//,6)(^i, ^^2) = d{^G,a){'"i) ■ d{^H,h){v2) 

Proof : This follows directly from equation (i) in Lemma 4.2. ■ 
Remark 4.4 : Prom corollary we get d(G,a)\i]{H,h) = d(^G,a) ■ d(H,b) 

Theorem 4.5 : Consider a bipartite system in W (g) in the state a. Then a — a2 
if and only if a is the density matrix of the graph (G, a) B [H, b), where (G, a) and {H, b) 
are the graphs having density matrices cxi and (72 respectively. 
Proof : If peirt : Given (G, a), {H, b) we want to prove 

a((G, a) B {H, b)) = ai(G, a) ® a2{H, b). 

Prom the definition of the modified tensor product we can write 

a((G, a) B {H, b)) = {Q[C{G, a) ® Crj{H, b) 

d(G,a)B{H,b) 

+C{G, a) ® Af{H, b) + Af{G, a) ® C{H, b) + f)(G, a) ® n{H, b)]} 
Using Lemma 2.10, Remark 2.11 and Remark 4.4 we get 

a{{G,a)B{H,b)) = ^- [Q{C{G,a) C7^{H,b)) 

d(G,a) • d(H,h) 

+Q{C{G, a) ® J\f{H, b)) + Q{U{G, a) 

0C(H,b)) + Q{n(G,a)0n(H,b))] (26) 
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We can calculate every term in (26) using (21) - (24) and substitute in (26) to get 

a{{G, a) B {H, b)) = a{G, a) ® a{H, b). 

Only if part : Given a = ai <^ a2, consider the graphs (G, a) and {H,b) for ai and (72 
respectively. Then the graph of a has the generalized Laplacian 

[L{G, a) + Ao{G, a)] ® [L{H, b) + Ao(i/, b)] 

= L(G, a) ® L{H, b) + L(G, a) ® Ao(G', a) + Ao(G, a) ® L{H, b) + Ao(G, a) ® Ao(i/, 6) 

Now it is straightforward to check that the graphs corresponding to each term are given 
by the corresponding terms in the definition of (G, a) □ (if, b) . ■ 
Remark 4.6 : Note that the proof of Theorem 4.5 does not depend in any way on the 
positivity or the hermiticity of the associated generalized Laplacians. Therefore we have 

Q{{G,a)B{H,b)) = Q{G,a)^Q{H,b) 

for any two graphs {G, a) and {H, b) 

Corollary 4.7 : The modified tensor product is associative and distributive with respect 
to the disjoint edge union +. 

Proof : Let (Gi, ai), (G2, 02) and (G3, 03) be any graphs. Using Theorem 4.5 and Remark 
4.6 , we can write 

g(((Gi,ai)H(G2,a2))H(G3,a3)) 

= Q((Gi,ai)H(G2,a2))®Q(G3,a3) 
= (Q(Gi,ai)0Q(G2,a2))®g(G3,a3) 
= (3(Gi,ai)0(g(G2,a2)(8)Q(G3,a3)) 

= g(Gi,ai)®g((G2,a2)H(G3,a3)) 

= Q((Gi,ai)H((G2,a2)H(G3,a3)) 

Therefore , 

((Gi,ai)H(G2,a2))H(G3,a3) = (Gi, ai) □ ((G2, 02) □ (G3, 03)) 

Similarly, using Lemma 2.10 and distributive property of the matrix tensor product 
we get 

g((Gi,ai)H((G2,a2)) + (G3,a3))) 
= g((Gi, ai) □ (G2, 02)) + ((Gi, ai) □ (G3, 03)) 
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Which gives 

{G,,a,)B{{G2,a2)) + {G,,a,)) = (d, ai) □ (G2, as) + (Gi, ai) □ (G3, ^3) 

Definition 4.8 : The cartesian product of two weighted graphs {G, a) and {H, b) is 
denoted {G, a)n{H, h) with weight function c defined as follows. 

V{G,a) X V{H,h). 

E{{G, a)n{H, b)) = {{{u, v), (x, u = x and {v, y} e E{H, 6), v ^ ?/, c({(u, y), (u, y)}) = 
c^u ■ ^({"^i y}) or V = y and {-u, G a), u ^ x, c{{{u, v), {x, v)}) = ■ a({n, x}). 

Where and arc the degrees of the vertices u € E{G, a) and v € E{H, b) respec- 
tively. It is straightforward to check that 

{G,a)U{H,b) = C{G,a)®Af{H,b)+Af{G,a)(^C{H,b) 

Which can be taken to be the definition of the cartesian product of graphs in terms 
of the operators C and M . We also note that 

{G,a)B{H,b) = {C{G,a)^Cr}{H,b)} + {{G,a)D{H,b)} + {n{G,a)^fl{H,b)} 

Note that the isolated vertices in (G, a) or {H,b) do not contribute to {G,a)D{H,b) 
as their degree is zero. 

Example(6) : Consider {G,a),{H,b) where V{G,a) = {1,2}, E{G, a) = {{1,2}} and 
V{H, b) = {1, 2, 3, 4}, E{H, b) = {{1, 2}, {2, 3}, {3, 4}} with weight functions a = 6 = 1, 
as shown in Figure 14a, 146. The modified tensor product of these graphs is given by 
Figures 15a, 15b, 15c, 15d, for each term in (25) and the resulting graph is as shown in 
Figure 15e. The corresponding density matrix of the graph {G, a) □ {H, b) is 

1 -1 -1 1 " 

-12-10 1-2 1 

0- 12-10 1-21 
-1 1 1 -1 

-1 1 1 -1 0' 

1- 2 1 0-12-10 
1-2 1 0-12-1 
1 -1 -1 1 _ 

which is the same as a{G, a) ® cr{H, b) . 
Corollary 4.9 : The density matrix of the modified tensor product of two graphs is 
separable. 



a{{G,a)B{H,b)) = ^ 
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Proof : Prom Theorem 4.5 we see that cr((G, a) □ {H, b)) is actually a product state. ■ 
CorollEiry 4.10 : cr = (7i (g) (T2 (8) • • • (8) (7jk for a /c-partite system if and only if the graph 
of a is the modified tensor product of the graphs of cti, • • • , (7^. 

Proof : Apply Theorem 4.2 successively to {ai (T2), (((Ji (8) 0-2) <8) (73) • • • and then use 
the associativity of the modified tensor product corollary 4.7. ■ 
Corollary 4.11 : A state cr of a /c-partite system is separable if and only if the graph 
(G,a) for a has the form 

{G,a) = U,Hj=,(Gi,a]). 

Proof : Let a be separable i.e. 

i i 
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By Algorithm 2.9 and Corollary 4.10 the graph of a has the form 

Now let the graph of a /c-partite state be 

(G,a) = U, {Gi,ai). 
Then by Lemma 2.10 , Remark 2.11 and the above corollary to Theorem 4.5 



a{G, a) = ^Wjcrj 



(1) 



.(2) 



cr- 



(fe) 



if the 



Corolary 4.11 says that Werner's definition [1] of a separable state in M'^^ (8)M'^2 (gjJR'Js , 
• • • ® M'^'= system, can be expressed using corresponding graphs. 
Lemma 4.12 : For any n = pq the density matrix a{Kn, a) is separable in R^(g 
weight function is constant > 0. 

Proof : The proof is same as that given for corollary 4.3 in [7], for simple graph. ■ 
Example (7) : Consider the graph (^4,0). The vertices of {K4,a) are denoted by 1, 

2, 3, 4, where weight function is constant, say , a = 3 > and has loops in vertices 1, 
2. We associate to these vertices the orthonormal basis {|1) = |1)|1), |2) = |1)|2), |3) = 
|2)|1), |4) = |2)|2)}. In terms of this basis <t{K4, a) can be written as 



12 
-3 
-3 
-3 



-3 
12 
-3 
-3 



-3 
-3 
9 
-3 



-3 
-3 
-3 
9 



1 

14 



4 
-1 
-1 
-1 



4 
4 
-1 
-1 



-1 -1 
-1 -1 



3 
-1 



-1 
3 



and from equation (15) we can write a^K^, a) as 
a{K4, a) -■ 



l[6P[|l)i=(|l) - |2))] + 6P[-^(|1) - |2))|1) 



+6P[-^(|1) - |2))|2)] + 6P[|2)i=(|i) - |2))] + 6P[-J=(|11) - |22))] 
6P[i=(|l2)-|21))] + 3P[|ll)] + 3P[|12)]} 

^(^4,a) = ip[|l)-L(|i)_|2))] + ip[-L(|i)_|2))|l)] 
+\p[^{\l)-\2m]+^-P[\2)^{\l)-\2))] 
+^{lp[i=(|ll) - |22))] + IPl^m) - |21))]} 
^^[|11)] + ^^[|12)]. 
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Each of the first four terms in the above expression is a projector on a product state, 
and also the last two terms are projectors, while the fifth and sixth terms give rise to the 
separable density matrix ^p[\ — ) | +)] + )], where |±) ± |2)) [7]. Thus 

a{K4^, a) , a constant, is separable in ® M^. 

Note that there exists a graph which is complete with a real weight function, which is 
entangled as the following graph shows in 




Figure 16 



Remark 4.13 : Separability of (T{Kn,a) with constant weight function > does not 
depend upon the labeling of V{Kn, a) provided every vertex has a loop or there are no 
loops. Given a graph, an isomorphism from (G, a) i — > (G, a) is called automorphism. 
Under composition of maps, the set of automorphisms of (G, a) form a group denoted 
Aut{G,a). If a{Kn,a) is separable, and if the Aut{Kn,a) = Sn, (G, a) = {Kn,a) is also 
separable. Note that Aut{Km a) — Sn provided all weights are equal and either every 
vertex has a loop or there are no loops. 

Lemma 4.14 : The complete graph {Kn, a) on n > 2 vertices corresponding to a separable 
state with weight function constant > is not a modified tensor product of two graphs. 
Proof : It is clear that, if n is prime then {Kn, a) is not a modified tensor product of 
graphs. We then assume that n is not a prime. Suppose that there exist graphs (G, b) 
and {H,c) respectively on p and s vertices such that {Kps,a) = (G, 6) □ {H,c) where b 
and c are constants . From the definition of the modified tensor product 

a({(ui, vi), (u2, V2)}) = b({ui, U2}) ■ c({vi,V2}), 

the degree sum is 

dicb)^ I] Kv^2b\E{G,b)\. 

ueV{G,b) ueV{G,b) weV{G,b) 

We know that d(^G,b) < b{p{p — 1)) and also di^H,c) — '2c\E(H, b)\ < cs{s — 1) and 

d{G,b) ■ d{H,c) < bcps{p - l)(s - 1) = bcps{ps -p- s + 1). (27) 
Now observe that V{{G, b) □ (if, c)) = ps and, 

d{G,b)B{H,c) = aps{ps - 1) (28) 
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because {G, h) □ (ff, c) — {Kps, a). 
We know that 

d(G,b)a{H,c) = d(G,b) ■ d(H,c)- (29) 

Substituting from (27) and (28) we see that (29) is satisfied only when p — 1 — s , i.e. 
n = ps = 1. ■ 

Lemma 4.12, Lemma 4.14 and Theorem 4.5 together imply that a complete graph 
{Kn, a) on n > 2 vertices with a = constant > is a separable state but not a product 
state. 

Definition 4.15 : Consider a graph (G, a), without loops, pertaining to a bipartite 
system of dimension pq . The partial transpose of (G, a), denoted {G^'^,a') , is a graph 
defined as V{G^'',a') = V{G,a) , {il,kj} e E{G^'',a') <^ {ij^kl} e E{G,a) and 
a'{{il,kj} = a{{ii, kl}. 

Lemma 4.16 : Consider a bipartite separable state (t(G', a) with the associated graph 
(G, a) without loops. Then A(G, a) = A(G'"^, a'), where {G^^ , a') is the partial transpose 
of (G, a). 

Proof : Let Q{G, a) be the Laplacian of a graph (G, a) with real weights without 
loops, on n vertices. Let D be any nxn real diagonal matrix in the standard orthonormal 
basis ^ = 1, 2, . . . , n, such that D and Tr{D) — 0. This means that there is 

at least one negative entry in the diagonal of D. Denote this element by — 6j. Let 
iV'o) = l^j) and 10) = Ej Xjkj) where, 

r iij^i 

\keR ilj^i 

Let Ix) = IV') + 10) = E:Li(1 + X.)k.)- Then 

{x\Q{G,a) + D\x) = (V'o|Q(G, 0)1^^0) + (V'o|Q(G,a)|0) + (0|g(G, 0)1^^0) + 

(0|Q(G,a)|0) + (VoPlV'o) + (V'oPI^) + {<t>\D\i^o) + (0P|0) 

Since l-i/^o) is (unnormalized) vector having all components equal unity, from equation (9) 
it follows that (V^o|(3(G, a)|^o) = 0. Also (V'ol^l^o) = Tr{D) = 0. We have 

{(t>\Q{G,a)\(t>) ^ k''{Q{G,a))u^ kHi 

(V^o|g(G,a)|0) = (0|Q(G,a)|^o) = O. 
Finally, the remaining terms in the above equation are given by 

(01^10) = h,k^ 
{MD\(t>) = hk^{4>Wo). 
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Thus 

{X\Q{G, a) + D\x) = k^ih + di) + 2kk 
So we can then always choose a positive k , such that 

{x\QiG,a)+D\x)<0. 

It then follows Q{G, a)+D^O. 

This expression is identical with that obtained in [2]. For any graph G on n — pq 
vertices 

Vl = UlWl, V2 = U1W2, ...,Vpg = UpWg, 

consider the degree condition A(G') = A{G^'^). Now 

{L{G)f^ = {A{G) - A(G^«)) + L{G^^). 

Let 

L> = A(G') - A(G'^^). 
Then D is an n x n real diagonal matrix with respect to the orthonormal basis 

\Vi) = \Ul) (g) \Wl), \Vpg) = \Up) (g) \Wg). 

Also 

Tr(D) = Tr(A(G)) - Tr(A(G^'')) = 0. 

We have two possible cases : D ^ or D — 0. If D 7^ 0, that is the degree condition is 
not satisfied (i.e.A(G) 7^ A{G^^)) we have seen that L{G) + D ^ 0. As a consequence, 
L(G'rs) + D ^ and then {L{G)f'' ^ 0. Hence p{G) is entangled. ■ 

Lemma 4.17 : A graph (G, a) for a bipartite state corresponds to a separable state 
if {ij, kl} {i k,j ^l) e E{G, a) =^ {il, kj} e E{G, a) and aij^ki = au^kj- 

Proof : Suppose Qij^ki — 0'ii,kj — (1,^ k, j ^ 1. The contribution of the corresponding 
two edges is 

a{p[^m-\ki))]+p[^{\u)-\km 

which is a separable state. Thus all such pairs contribute separable states. Any other 
edge {ij, kl} with i — k or j — I has the contribution aij^kiP[\i) ® (71 (li) ~ 10))] which is 
separable. Loops contribute the product states P[|m)]. ■ 
The reverse implication is not true in general. The counter-example is the graph 
(Figure 12a) in example (4) which is separable. 
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5 Graph Operators 



A graph operation is a map that takes a graph to another graph [17]. We deal with four 
cases namely deleting and adding an edge and deleting and adding a vertex. 
Deleting an edge {{vi, Vj}, ay^y.)] from a graph {G, a) results in a graph 

def 

{G, a) — {{vi, Vj}, av^vj) — a), E{G, a) \ {vi, Vj}) with av^vJ — 0. Note the possibility 

Vi — Vj corresponding to the edge being a loop. Addition of an edge {{ui,Vj},aij) maps 

def 

{G, a) to the graph {G, a) + {{vi, Vj}, Uij) — [V{G, a), E{G, a) U {vi, Vj}] with a^.„^. = aij. 

def 

Deletion of a vertex Vi maps {G, a) to {G, a) — {vi} — [V{G, a) \ {vi}, E{G, a) \ Ei] where 
Ei is the set of all edges incident to Vi (including the loop on Vi) with the weight function 

def 

zero for the edges in E^. Adding a vertex Vi to {G, a) maps {G, a) to {G, a) + {vi} — 
{V{G,a)U{vi},E(G,a)). 

A very important point is that, in general, the set of graphs having density matrix is 
not closed under these operations. Addition of an edge with positive weight and deletion 
of an edge with negative weight preserves the positivity of the generalized Laplacian 
resulting in the graph having density matrix. However, addition (deletion) of an edge with 
negative (positive) weight may lead to a graph which does not have density matrix. In the 
next section, we give a method for addition and deletion of an edge which preserves the 
positivity of the generalized Laplacian. Deletion and addition of vertices always preserves 
the positivity of the generalized Laplacian. 

Let i3(7Y") be the space of all bounded linear operators on Ti". A linear map A : 
B{T-C"') B{T-C"^) is said to be hcrmiticity preserving if for every hermitian operator 
O G i3(7i"),A(0) is an hermitian operator in i3(7i™). A hcrmiticity preserving map 
A : B{W^) B{W^) is said to be positive if for any positive operator O E B{W^),A(0) 
is a positive operator in i3(7Y™). A positive map A : BiTi^) — > B{H^"') is said to be 
completely positive if for each positive integer k, (A ® 1^2) : i3('H" Ti^) B{H'^ H'^) 
is again a positive map. A completely positive map A : B(7i") — > B{VJ^) is said to be 
trace preserving if Tr{h.{0)) — Tr{0), for all O e B{TiP-). A quantum operation is a 
trace preserving completely positive map (for short, TPCP) [9, 5]. In standard Quantum 
Mechanics, any physical transformation of a quantum mechanical system is described by 
a quantum operation [6]. We are going to use the following result: 
(Kraus representation Theorem) [10] : Given a quantum operation A : i3(7i"') 
B(TC"^), there exist mx n matrices Ai, such that A(p) = ^ A^pAl, where p is any density 

i 

matrix acting on and X^^j^A = Im (The converse is true). The matrices AiS are 

i 

called Kraus operators. 

A projective measurement M. — {Pi, i — 1,2, ■ ■ ■ , n}, on a quantum mechanical sys- 
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tern S whose state is p, consists of pairwise orthogonal projectors Pi : Hs ^ Hs, such 

n 

that Pi — Idim{Hs)- The i-th outcome of the measurement occurs with probabihty 

i=l 

Tr(Pip) and the post-measurement state of S is ^^f^- Whenever the i-th outcome of the 
measurement occurs, we say that Pi chcks. Last two paragraphs apply to complex Hilbert 
space and so also to real Hilbert space. 

5.1 Deletion and addition of an edge for a weighted graph with 
all weights > 

Here we describe how to delete or add an edge by means of TPCP. Our method of 
deleting an edge from a weighted graph with all positive weights is a simple generaliza- 
tion of the method in [7]. Let {G,a) be a graph on n vertices • ■ • and m edges 
{vijVj^} ■ ■ ■ {vi^Vj^}, ik 7^ jk,k = 1, - ■ ■ ,m and s loops {vi^Vi^} ■ ■ ■ {vi^ViJ. Our purpose 
is to delete the edge {vii^Vjj^}, ik ^ jk- Then we have 

m ^ s 

a{G,a) = {J^2ai,j,P[— =(|vi,) - It;^,))] +^ai^i^P[\vi^)]} 

and 



a{{G,a) - {vi^VjJ) - 



d{G,a) - '^^ikik 



m ^ s 

1] '^'^imP[-^i\'"k) - 1%))] + ai^kPiM] 



> . 



A measurement in the basis Ai — {^{\vi^,)^\vj^)), \vi) : i ^ ik,jk a-nd i — 1,2, ■ ■ ■ ,n} 
is performed on the system prepared in the state a{G,a). The probability that P+ = 
P[^(K) + 1%))] clicks is 



Tr[PMG,o,)] = Y,{^i\PMG,o,)\vi 



1=1 



^ m s 

77j { / ^ '^iejA^ikH ~ ^ikk "I" ^jkii ~ ^jkji\ + / ^ '^itH {^Wk + ^itjk) } 



(30) 



The state after the measurement is P[^{\vi^) + \vj^)]. Let f/^ and f/^ be nxn unitary 
matrices such that = -^{\vit) -\vj^)) for £ = 1, • • • ,k-l,k + l, - ■ ■ ,m 
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and U^t[^{hk) + \vjk))] = K),t =!,■■■ ,s. Now, with probability 2ai^jJ{d(^G,a)-2ai^jJ 
we apply f/^ on P[-^(\vi^) + \vjj] for each i = 1, ■ ■ ■ , k — 1, k + 1, ■ ■ ■ ,m .and with 
probability ai^iJ{d(G,a) - ^ai^jj we apply f/^ on P[-^{\vi^) + \vj^^)] for each t = 1, - ■■ ,s. 
Finally we obtain o-{{G, a) — jj,}) with probability given by (30). The probability that 
P[^iK) - K))] clicks is 



2d 



i = i 



{G,a 



(31) 



the state after measurement is P[-^{\vii^) — I'fjk))]- and Uj^^ nxn unitary matrices 

such that ^ 

for £ = 1, • • • , A; — 1, A; + 1, • • • , m and 



Ukt^iK) - 1%)) = K) 

fort = 1, -,5. With probability 2ai,j,/(rf(G',a)-2aifcjJ we apply L^^^ on P[^(|t'iJ-|t'jJ)] 

for each i = ,k — l,k + l,-- - ,m and with probability ai^ij {d(^G,a) — '^(^ikjk) *c apply 

f/^7i on P[^{\vi^^) - \vji^))] for each t = 1, 2, • • • , s. Finally we obtain a{{G, a) - {vi^vj^}) 
with probability given by (31). 

The probability that -P[|'yi)] where i ik,jk and i — 1, - ■ ■ ,n clicks is 



(G,a) 



i=l 



t=i 



(32) 



and the state after measurement is P[|vj)]. Let Ua and Uu he n x n unitary matrices 
such that Uii[\vi)] = ^(It'i^) - \vj^)] for i = 1, ■ ■ ■ k - 1, k + 1, ■ ■ ■ ,m and Uit[\vi)] = \vi^) 
for t = 1, ■ ■ ■ , s. With probability 2ai^jJ{di^G,a) — '^cudk) apply Uu on for each 

£ = I,-- - ,k — \,k + 1, ■ ■ ■ ,m and with probability ai^ij {d(G,a) ~ '^C'ikjk) apply Uu on 
P[|vi)] for each t = 1, ■ ■ ■ , s. 

We obtain a{{G,a) — {vi^.,Vj^}) with probability given by (32). This completes the 
process. 
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The set of Kraus operators that reahzes the TPCP for deleting the edge {vi^,Vji^} is 
then 



{ 



2a,,, U^^P[^{\v,^) + \v,^))];i =l,...,k-l,k + l,---,m} 
a{G,a) - ^aikjk v2 



"nn __L„. 1 



a{G,a) - ^aidk V 2 



lj| / 2ai^j, 1 



UuPi^iK) - K))] ■.£ = l,---,k-l,k + l,---m} 



U{ 



u{, 



u{ 



2a,; 



-UitP[\ui)] : i = 1, • • • , n, i 7^ ife, jit; £ = 1, • • • , A; - 1, A; + 1, • • • , m} 



UitP[\vi)] : i ^ 1,- ■ ■ ,n,i ik,jk;t = 1,- ■ ■ ,s} 



The set of Kraus operators that reahzes TPCP for adding back edge {vif.,Vji^} to (G, a) 



a(G,a) + ^aikjk V 2 



U 



U 



U 



U 



U 



a{G,a) + mujk V 2 

;r^5|— - 1%))] : ^ = 1, 2, • • • , m 

"(G,a) + m^jk V 2 



-VuP[\vi ^] : i = 1, 2, • • • , n, i 7^ ijk, jjfc, £ = 1, 2, • • • , m 



d{G,a) + 2a, 



-VitP[\vi)] : i = 1,2, • • • 7^ ijfc, = 1,2, • • • ,s 
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where VVt, Vl«, , V^, Vu are n x n unitary matrix defined as follows: 



1 



1 



'^ie) - \voi)), for ^ = 1,2, • • • ,m 



^fet;^(l^iJ + \'^h)) = W)^ for i = 1, • • • , s, 
1 1 



- 1%)) = l^it))' for t = I,-- - ,s 



1 



V'i^ki) = -^(ki,) - ki,)), for i 1,2, • • • 7^ ifc,jfc,£ = 1, • • • ,m 

^itki) = kit), for i = 1, • • • 7^ 4, = 1, • • • ,s. 

For deleting a loop {t'i^, , Vj^, } a measurement in the basis = 1, ■ ■ ■ ,n} is per- 

formed on the system prepared in the state (7(G,a). Then the probability that -P[ki)] 
clicks for i = 1, • • • n is 



^ m s 



t = 1 
t^t' 



(33) 



The state after the measurement is P[|vi)]. Let Uu be n x n unitary matrices such that 
U-d[\vi)] = -j^i\vi,) - \vj^)). For i = I,-- - ,m and Uit[\v,i}] = \vi^), for t = 1, ■ 



l,t' + 1, 



t' - 



,s. With probability 2ai^jJ (d(^Q^a) ~ '^iti,iti) wc apply on -P[|fi)] for each 



£ = 1, ■ ■ ■ ,m and with probability cutitl {d{^Q^a) — cUtiiti) we apply Uu on ^[It'i)] for each 
t = 1, ■ ■ ■ ,t' — l,t' + 1, ■ ■ ■ s. We obtain cr((G, a) — {^^^,^1^,}) with probability given by 
(33). 

The set of Kraus operators that realizes the TPCP for deleting the loop {vi , , Vi , } is 



\G,a) 



-UuP[\vi)] i = l,---,m, £ = 1, 



, m 



U 



2a, 



d 



(G,a) 



-UitP[\vi)] i = l,---,m, t = l,--- ,t' -l,t' + l,---s 



The set of Kraus operators that realizes the TPCP for adding the loop {vi^,Vi^,} 

-VieP[\vi)] : i = 1, • • • , n, £ = 1, • • • , m 



2a. 



d(G,a) + 
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^Sri^l^^VitP[\v,)]:i^l,---,n,t^l,---,s] 
where Va, are n x n unitary matrices define as follows : 

Vu\vi) = -^{\vii,) - \vj^)), for £ = 1, • • • ,m,i = 1, • • • ,n 
Vit\vi) ^ \vi^) , for i = 1, • • • = 1, • • • ,n. 

5.2 Deletion and addition of an edge with real weight, which 
preserves the positivity of the generalized Laplacian 

Let (G, a) be a graph with real weights on its edges not necessarily positive. We are 
basically concerned here with the deletion of {vi,Vj} with a^.y. > and the addition of 
{vi, Vj} with av^y^ < 0, because in other cases the positivity of the Laplacian is preserved. 
We define the sets 

= {{vi, Vj} e E{G, a), a,,y. > 0}, (34) 
E- = {{vi, Vj} e E{G, a), a„,„, < 0} (35) 

SindE^E+UE-. 

We define a graph operator S as 

E[E] = EU {{vi,Vi}, {vj,Vj} : a^.^. = a^.y. = 2|a^.^J and {vi,Vj} e E~} (36) 

Suppose we wish to delete a positive weighted edge {vii^,Vji^} e E'^ then we define the 
resulting graph as 

EC{{G,a) - {vi„Vj^}) 

where the graph operator jC is defined in (20b). 

For adding a negative weighted edge between Vi and Vj. i ^ j, we act on E{G,a) by 
the appropriate element of the set of operators {&ij}, i, j = I, ■ ■ ■ ,n,i ^ j defined as 

Eij [E]^ E\J {{vi, Vj}, {vi, Vi}, {vj, Vj} : a„,^^. < 0, a^,„, = 2|a„,„^. | = a^^.„J (37) 

To obtain the set of the corresponding TPCP operators we decompose the resulting graph, 
{G',a') given by EC{{G,a) - {vi^^,VjJ){ay^^y^^ > 0) (eq. (36) ) or by Gij ((G,a) + 
{vi,Vj}){ay^y^ < 0) (eq. (37)) or by {{G,a) - {ui^VjJ){ay^^y^^ < 0) or by (G, a) + 
{viVj}){ay^y. > 0) into spanning subgraphs determined by the sets E'^ and E" and 
treat the spanning subgraph corresponding to E~ replace the weights a^.^^. of edges 
{vi,Vj} G E~ by —tty^yj, so that both the spanning subgraphs have only positive weights. 
For getting the Kraus operators we go through the following steps. 

(a) First we determine the degree sums for the resulting graphs {G',a') in four cases. 
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(i) Deletion of a positive weighted edge {uii^jVj,^} 

d{G',a') = d(G,a) " 2a^i^^^ + 2 \ay^y. \ (38) 

i {ui,Vj}eE- 

(ii) Addition of a positive weighted edge {vi,Vj}. 

d(G',a') = d(G,a) + 2a^.^^ (39) 

(iii) Deletion of a negative weighted edge {vi^ , Vj^ } 

d(G',a') = d(^G,a) " 2a^i^«,.^ (40) 

(iv) Addition of a negative weighted edge {vi,Vj} 

d{G',a') = d(^G,a) + '^0-Vi,Vj + ^\o-Vi,Vj \ (41) 

(b) We construct the Kraus operators separately for and G~ for deleting the same 
edge {f j, Vj} from U {vj, Vj} or adding the edge {vi, Vj} to G^, using the method given 
in Section 5.1. However, the probability of applying various unitary operator [/^^ and 
^kli ^^'^ determinate using di^G',a') as in step (a) above. 

(c) Let {A.j} and {Bi} denote the sets of Krous operators for the graph operations on 
and G^ as described in (b). Then 

a{G\ a!) = ^^'^(^' «) - E ^^■^(^' (42) 

i 3 

and 

Y.A^^-Y.B]B,^I (43) 

« 3 

which can be justified by construction. 

We comment here that it is possible to modify the graph, after deleting a positive 
edge or adding a negative edge, which can preserve positivity in different ways, leading to 
different sets of Kraus operators. The basic idea is to add new loops. In our method we 
try to minimize the addition of loops. Further, in our method we cannot reverse the graph 
operation for deleting a positive edge or adding a negative edge. But this is not a problem 
since the quantum operations given by super operators are, in general, irreversible. 
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5.3 Deleting Vertices 

In order to delete a vertex Vi from a graph (G,a), 

(i) Delete edges, including loops, on Vi, one by one, by successively applying the proce- 
dure in 5.2. The resulting graph (G", a') has density matrix with i-th row and i-th 
column containing all zeroes. 

(ii) We now perform, on (t(G", a'), the projective measurement M = {/„— P[|f P[|wj)]}. 
Since /'[jfj)] is the matrix with all elements zero except the i-th diagonal clement, 
while o"(G", a') as all zeros in i-th row and column,the probability that -P[|fi)] clicks 
= Tr{aP[\vi)]) = 0. Thus when M is performed on a{G', a')', In-P[\vi)] clicks with 
probability one and the state after measurement is a{G', a') — {vi}) and is the same 
as a{G',a') without i-th row and i-th column. 

Adding a vertex : Let (G, a) be a graph on n vertices f i, ■ ■ ■ ,Vn and m edges {fi^fj^}, k = 
1, • • • ,m,ik 7^ jk and s loops {v^^Vi^}, t = 1, - ■ ■ , s; < ik,jk, H < n. Consider the following 
density operator 

P= {^^l^^^n\u^)])®{a{G,a)^ 

where \u2)} form an orthonormal basis of C^. We associate vertices Ui,i = 1,2 to 

the state \ui). Consider the graph H — ^2}, {{^i, ^i}, {^25 ^2}}) with associated 

2 

weights 6 > 0. It is easy to check that cr(if, &) = | ^ biiP\\Ui)]- Also observe that 

i=l 

p - a{{H, h) □ (G, a)) = a(H, b) ® a{G, a). 

Thus {H, b)E\{G, a) is the graph on 2n vertices labeled by UiVi. ■ ■ ■ , UiVn, U2V1, ■ ■ • , U2Vn 
and with 2m edges and 2s loops (see Section 4.2) {uiVi^.UiVj^} ■ ■ ■ {uiVi^.UiVj^^} 
{u2Vi^,U2Vj^} ■ ■ ■ {u2Vi^,U2Vj^} and loops {uiVi^,UiVi^},{u2Vi^,U2Vi^},t = l,---,s. So 
{H, h) □ (G, a) = (i^i, ai) W {H2, 02) where 

{Hi, ai) = {{uiVi ■ ■ ■ UiVn}, {{uiVi^,UiVj^} ■ ■ ■ {uiVi^,UiVj}}) 

{H2,a2) = i{u2Vi ■ ■ ■ U2Vn}, {{u2Vi^,U2Vj^} ■ ■ ■ {u2Vi^,U2Vj^}}) 

We first delete all edges and loops of {H, b) □ (G, a) which are incident to the vertex 
U2V1 G V{H2, 02) as in Section 5.2. Now we perform the following projective measurement 
on a{{H, b) □ (G,a)) 

n n 

M = {hn -Y,P[\u2V,)iY,P[M\} 

i=2 i=2 
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The probability that — X] ^[^^2^1)] is one and the state after the measurement is 

i=2 

a((Hi,ai) + {u2Vi}). 

Example (8) : Consider the graph as given in the Figure (16) , we want to delete the 
edge {1,2} with positive weight by means of TPCP. calculate the Kraus operators for 

and G~ as in section 5.2 ,where Ai for i — 1, • • - 24 and Bi for G~ ,i — 1, ,4 and 

substitute in the the equation 

a{G', a') = J2 A(^{G, a)4 - Yl ^M^, a)B] 

i 3 

where 

1 -1-1 1 " 

-1 3 -1 -1 

-1 -13-1' 

1 -1-1 1 

we get 

2 0-11" 

2-1-1 

-1 -13-1' 
1-1-13 

and 

24 4 

1=1 j=i 




Figure 16 



6 Representation of a general hermitian operator by 
a graph 

In this section, we generalize sections 2 - 4, to quantum states in a complex Hilbert space, 
that is, to the density matrices with complex off-diagonal elements. We have also given 
rules to associate a graph to a general hermitian operator. We believe that any further 
advance in the theory reported in this paper will prominently involve graph operators and 
graphs associated with operators. 



a{G,a) = l 



(^{G,a) = — 
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6.1 Representation of a general density matrix with complex off 
diagonal elements 

Consider a. n x n density matrix with complex off-diagonal elements. We associate with 
this density matrix an oriented graph {G,a) on n vertices, m edges and s loops with 
weight function 

a : V{G) X V{G) C. 

The weight function a has the following properties: 

(i) a{{u,v}) 7^ if ?'} e E{G,a) and otherwise. 

(ii) a{{u,v}) = a*{{v,u}) 

wc write a({-u, f}) = |a({-u, f })| e*''^"'', 0^,^, = 0. 

Note that, when (pij = InJ = 0, 1, • • • , i.e. a{{u, v}) is real, positive when / is even 
and real negative when / is odd. 

The degree of vertex v is given by 

d{G,a){v) = o?„= J2 + (44) 

u&V{G,a),u^v 

d(G,a) = ^ dy 
v&V{G,a) 

The adjacency matrix M(G', a) of a complex weighted graph with n vertices is a n x n 
matrix whose rows and columns are indexed by vertices in V{G, a) . 

M„„ = a{{u,v}) = arUv^u}) = (M„„)*. 

The degree matrix A(G', a) of the complex weighted graph is a n x n real diagonal 
matrix, whose rows and columns are labeled by vertices in V{G, a) and whose diagonal 
elements are the degrees of the corresponding vertices. 

A(G, a) = diagldy] v e V{G, a)] 

where (Z„ is given by equation (44). 

The loop matrix Ao(G', a) of a graph (G, a) is a n x n real diagonal matrix with 
diagonal elements equal to the weights of the loops on the corresponding vertices 

The generalized Laplacian of a graph (G, a), which includes loops, is 

g(G,a) = A(G,a) + M(G,a) - Ao(G',a) (45) 
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Note that Q{G, a) is hcrmitian matrix . If the generahzed Laplacian Q{G, a) is positive 
semidefinite, we can define the density matrix of the corresponding graph {G, a) as 

1 



a{G, a) 



(G,a) 



Q{G,a) 



(46) 



where Tr(a(G,a)) = 1. 

For any n x n density matrix a with complex off diagonal elements we can obtain the 
corresponding graph as follows: 
Algorithm 6.1 : 

(i) Label the n vertices of the graph by the kets from the standard orthonormal basis. 

(ii) For every nonzero ijth element with j > i given by a{{i,j}) draw an edge between 
vertices labeled \vi) and \vj), with weight a{{i,j}). 

(iii) Ensure that dy^ = an by adding loop of appropriate weight to Vi if necessary. 
Example (9) : (1) 

11 



(a) 



1 



(b) 

Figure 17 




where |y, +) 
(2) 



where |y, 
(3) 



PU+)] 



PU -)] = 2 



" 1 -i " 


1 


1 e-'^/2 - 




i 1 


~ 2 


g-m/2 I 




\e correspondin 


g graph is as shown 


1 i 


1 






-i 1 


~ 2 







;^(|1) — i\2)) and the corresponding graph is as shown in Figure 17b 



P[\y,+)\y,+)] 



1 -I ~i -1 

i 1 1 -i 

i I I ~i 

-1 i i 1 
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The corresponding graph is as shown in Figure 17c. 

Note that Remark 2.1 is vahd also for complex weighted graphs. 

Remark 6.2 : Theorem 2.3 applies to complex weighted graphs with equation (11) 
changed to 



m 



5^<2 + 2j]|a,,,J^ = 4G,a) (47) 

i=l fe=l 

also Lemma 2.4 applies to complex weighted graphs. 

Definition 6.3 : A graph [H, h) is said to be a factor of graph (G, a) if V{H, h) — V{G, a) 
and there exists a graph (//', h') such that V{H', h') = V{G, a) and M{G, a) = M{H, h) + 
M[H',h'). Thus a factor is only a spanning subgraph. Note that 



h^^^. if {vi.Vj] e E{H,h) 
b' if {v„Vj}eEiH',b') 



Now let {G, a) be a graph on n vertices Wi, ■ ■ ■ , w„ having m edges 
{vh:Vj^}, • • • , {vi^,Vj^} and s loops {v^^, v^J • • • {vj,, v^J where 1 < iiji, • • • , < 
n, 1 < iii2 ■ ■ - is <n. 

Let {Hi^ji^, ttii^j^) be the factor of (G, a) such that 

[Mm « ^1 -/ "^^^fc " = 't^^jk or a*^^.^if u^jk,w^ik . s 

FA/H- ^ A] _/ WiuJ if = = w or u^jk^w 
m-t^ikh,ai,jj\u,nj - <^ Q Otherwise ^ ' 

Let OitiJ be a factor of {G, a) such that 

\A/r(u W —\\(u w —f when u = if = w , , 

[M(iiijij, aj^ijju^tu = [A(//j^ij, ai(iJJ„,t„ = < ^ otherwise ^ 

Theorem 6.4 : The density matrix of a graph (G, a) as defined above with factors given 
by equation (48), (49) and (50) can be decomposed as 

^ m 1 ^ 

a{G,a) = V2|a({ife,jfe})|c7(i/'i^j,,ai,jJ + V ai^i,(T(i/i,i,, a^.^J (51) 

«{G,a) ^ a(G,a) ^ 



or 



m 1 * 

^(G,a) = J]2|a({^fe,j.})|P[-=(|^.J -e^'^'.^H^..))] + E "M^^'Ot'.)] (52) 
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Where (f)i^,j^, = vr for any edge {ik,jk} with real positive weight and (pi^j^. — for any real 
negative weight. 

Proof : From equation (48), (49), (50) and Remark 6.2, the density matrix 



is a pure state. Also, 



1 



di. 



"ItH 



is a pure state. Now 



A(G, a) = ^ A(i/i,j,, ai,jj + ^ Ao(i/i,j,, a^^ij 



k=l 



fe=l t=l 

Therefore , from eq. (46) 

1 



a{G,a) = 



fc=i 



k=l 



+ 



d 



(G,a) 



t=l 



d, 



^ JTt 1 ^ 



In terms of the standard basis, the uw-th element of matrices a{Hi^j^,ai^j^) and 
are given by {vu\(T{Hi^j^, , Ui^jJlv^) and {vu\cr{Hi^i^ai^i^\vw) respectively. In 



this basis 



<^iHikjk,aikjk) = ^[;7f (I^h) - (^"^''''K))] 



Therefore equation (51) becomes 



aiCa) = ^f]2|a(fejfc})|P[^(|^,,) -e^^^^^H^)) + 7^ E 

(52) 

Where (pi^j^ = tt for any edge {ik,jk} with real positive weight and ^j^j^ = for any real 
negative weight. 
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Example (10) : 

(1) For a graph given in Figure 17b, the density matrix is 



a{G,a) - 

(2) For a graph given in Figure 17c, the density matrix is 



a(G,a) = ^{2P[i=(|ll)-e-/^|12))] + 2P[i=(|ll)-e-^-/^|21))] 

+2P[i=(|ll) - |22))] + 2P[-^(|12) + |21))] 

+2P[i=(|12) - e-/^|22))] + 2P[^(|21) - e-/^|22))] 
-2P[|11)] - 2P[|22)] - 2P[|12)] - 2P[|21)]} 





1 


—i 


—i 


-1 


1 


i 


1 


1 


—i 


4 


i 


1 


1 






-1 


i 


i 


1 



Example (11) Consider the state 

1 



'^ = ^P[\y.+)\y.+)] + ^P[\y.+ 



where |y, +) = ^(|1) + 1\2)) and |^) = ^(|1) + iV2\2)) 



1 

36 



1 

36 



7 -(3 + 4\/2)i -7i -(3 + 4^) 

(3 + 4v^)i 11 3 + 4^2 -Hi 

7i 3 + 4^2 7 -(3 + 4V2)i 

-(3 + 4^2) Hi (3 + 4V2)i 11 

7 (3 + 4V2)e-^'^/2 7e-^^/2 -(3 + 4^2) 

(3 + 4 V2)e^'^/2 11 7 116-^^^/2 

76^^^/^ 3 + 4^2 7 (3 + 4 v^) 6-^^/2 

-(3 + 4^2) lle^'^/^ (3 + 4V2)e^^/2 



The corresponding graph is as shown in Figure 18 , 
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-(6+8K2) ,^ .J 

-(3+4|/2)i - 




(6+81/2") 



21 -(3+4|/^)i 22 



-(6+8K2) 



Figure 18 



and using the equation (45) to get the matrix from graph in Figure 18, 

a{G,a) = ^{2(3 + 4V2)P[-i=(|ll)-e-/2|12))] 

+2 X rP[-^(|ll) - e-/^|21))] + (3 + 4x/2)P[-L(|ii) _ |22))] 

+2(3 + 4V2)P[-i=(|12) + |21))] + 2 X 11^[;^(|12) - e-/^|22))] 

+2(3 + 4V2)P[^(|21) - e-''^/'|22))] - (6 + 8v^)P[|ll)] - (6 + 8V2)P[|22)] 
v2 

-(6 + 8V2)P[|12)] - (6 + 8y2)P[|21)]}. 
We can check that 



(t(G, a) 



1 

36 



7 -(3 + 4v^)i -li -(3 + 4x/2) 

(3 + 4^/2)i 11 3 + 4^/2 -Hi 

7i 3 + 4^2 7 -(3 + 4v/2)i 

-(3 + 4^2) Hi (3 + 4V2)i 11 



We can also check that this state is not pure by applying Remark 6.2 on the graph. 
6.2 Separability 

Remark 6.6: The definition of the tensor product (G, a)® (if, h) of two complex weighted 
graphs (G, a) and {H^h) is the same as given before. However note that {^1,^2} £ 
E{G, a), {wi, W2} e E{H, b) implies 



c{{{vi,Wi),{v2,W2)}) = a{{vi,V2})b{{wi,W2}) 



and 



c{{vi, W2), {v2, wi)}) = a{{vi, V2})b{{w2, wi}) = a{{vi, V2})b*{{wi, W2}). 
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Remark 6.7 : Equations (18a) and (18c) are valid for the tensor product of complex 
weighted graphs. Also, Q{{G,a) (g) {H,b)) ^ Q{G,a) ®Q{H,b). Equation (18b) holds 
good only for graphs without loops, for graphs with only loops or when one factor has no 
loops and other factor has only loops. For such graphs equation (18b) immediately gives 

d{G,a)®{H,b){v,w) = d^G,a){v) ■ d^H,b)(w). 

6.2.1 Modified tensor product 

The modified tensor product of two complex weighted graphs requires the operator J\f to 
be redefined in the following way. We replace the equation (20a) by 

= ^ \ai{vi,Vk})\ + a{{vi,Vi}) (53) 

Vk&V{G,a) 

The definitions of the operators r],/! and Q remain the same. Equations (21) to (24) are 
satisfied by these operators on the complex weighted graphs. We further have 

(i) M{^/C{G,a)) = A{G,a) - Ao{G,a) 
A{XCiG, a)) = AiG, a) - Ao(G, a) 
Ao{AfC{G,a)) = A{G,a)-Ao{G,a) 
QiAfCiG, a)) = A{G, a) - Ao(G, a) 

The modified tensor product of two complex weighted graphs {G, a) and {H, b) with p 
and q{> p) vertices respectively is 

(G, c) = (G, a) □ {H, b) = £(G, a) ® C{H, b) + C{G, a) ® Af{H, b) + Af{G, a) C{H, b) 

+{n(G, a) (g) n(H, b) U 2XjC(G, a) J\fCr]{H, b)} (55) 

The weight function c of (G, a) □ (if, h) is obtained via the definition of tensor product 
and the disjoint edge union. 

Lemma 6.8 : A((G, a) □ {H, b)) = A(G, a) ® A{H, b). 

Proof : Since Lemma 2.10 applies to disjoint edge union of complex weighted graphs, 
A((G, a) □ {H, b)) = A{C{G, a)®C{H, b)) + A{C{G, a) ^M{H, b)) + A{Af{G, a)®C{H, b)) 

+A((](G, a) ® Q{H, b)) + A(2AA£(G, a) ® MCri{H, b)) 

The last two terms are justified because the graphs involved are real weighted graphs. 
Using Remark 6.7 we get 
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A{{G, a)B{H, b)) = A{C{G, a))®A(£(iJ, 6))+A(£(G', a))0A{f^{H, b))+A{f^{G, a))®A{C{H, b)) 

+A{n{G, a)) (g) A{n{H, b)) + 2A{AfC{G, a)) (g) A{AfCr}{H, b)) 
Using equations (21) to (24) and (54), we get, after some simplification. 



A((G, a) B {H, b)) = A((G, a)) A{{H, b)) 

Corollary 6.9 : d^G,a)BiH,b)iv,w) = d^G,a){v) ■ d^H,b){w) 
and 

d{G,a)B{H,b) = d(^G,a) ' (^(if,6) 

Proof : The first result follows directly from Lemma 6.8. For the second note that 

Tr{A{{G, a) □ {H, 6))) = Tr(A(G, a) ® A{H, b)) = Tr{A{G, a)) ■ Tr{A{H, b)) 

where Tr denotes the trace. ■ 
Theorem 6.10 : Consider a bipartite syatem in C in the state a. Then 
cr = (Ti ® o"2 if and only if a is the density matrix of the graph (G, a) □ {H, b) where (G, a) 
and (if, b) are the graphs having density matrices ai and (T2 respectively. 
Proof : If part : Given (G, a), (if, 6) we want to prove 

(7((G, a) □ (ii, b)) = (7i(G, a) ® a2{H, b). 

From the definition of the modified tensor product we can write 

(7((G, a) □ {H, b)) = ^ {Q[C{G, a) C(H, b) 

a{G,a)B(H,b) 

+C{G, a) ^Af{H, b) +Af{G, a)^C{H, b) + {n{G, a)^n{H, b))U2AfjC{G, a) ^AfCr]{H, b)]} 
Using Remark 6.5 and corollary 6.9 we get 

a((G, a) B (H, b)) = ^--i^[Q(£(G, a) ® JC(H, b)) + Q{C{G, a) ® U{H, b)) 
+Q{M{G, a) (g) C{H, b)) + Q{yt{G, a) ® n{H, b) U 2MC{G, a) ® NCr]{H, b))] (56) 

We can calculate every term in (56) using (21) to (24) and (54) and substitute in (56) 
to get 

(7((G, a) B {H, b)) = a{G, a) (0 a{H, b). 
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Only if part : Given o" = o"i ® (T2 consider the graphs (G, a) and {H, h) for cti and 
(72 respectively. Then the graph of cr has the generalized Laplacian 



Q{G, a) ® Q{H, h) = (A(G', a) + M{G, a) - Ao{G, a)) (g) {A{H, h) + M{H, h) - Ao{H, b)) 

= A{G, a) A{H, h) + A(G', a) ® {M{H, h) - Ao{H, b)) + {M{G, a) - Ao{G, a)) ® A{H, b)+ 

{M{G, a) - Ao(G, a)) (g) {M{H, b) - Ao{H, b)) (57) 

Using equation (21) to (24) and (54) we see that RHS of equation (57) is the generalized 
Laplacian for {G, a) □ (if, b) ■ 
Remark 6.13 : The proof that the modified tensor product is associative and dis- 
tributive with respect to the disjoint edge union is the same as that for the case of real 
weighted graphs (corollary 4.7). 

Remark 6.14 : The definition of the cartesian product of graphs is the same as given in 
definition 4.8. 

Remeirk 6.15 : Corollaries 4.9 and 4.10 apply to complex weighted graphs without any 
change. 

6.3 Convex combination of density matrices 

Consider two graphs (Gi,ai) and (^2,^2) each on the same n vertices, having a{Gi,ai) 
and cr{G2, 02) as their density matrices respectively, where Oi and 02 are complex weight 
functions. Let (G, a) be the graph of the density matrix a{G, a) which is a convex com- 
bination of cr{Gi, Oi) and cr{G2, 02), 

a{G, a) = Xa{Gi, oi) + (1 - X)(t{G2, 02), < A < 1. 
It is straightforward, using the definitions of the operators J\f, L and 77, to varify that 

(G, a) = [AAr(Gi, ai) U (1 - A)A/'(G2, aa)] U [A£(Gi, ai) U (1 - A)£(G2, 02)] 

UriJC[XjC(Gi, ai) U (1 - A)£(Gi, 02)] (58) 
We can apply this equation to any convex combination of density matrices. Let 

a{G, a) = ^pia{Gi, a^), = 1 

i i 

Then, 

(G, a) = [U,p,Ar(G„ ai)] U [U,p,£(G,, a,)] 

Ur)C[UiPiC{Gi, ai)] (59) 
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where a and {flj} are complex weight functions, a{{vi, Vk}) — diiivi, Vk}) and a{{vi, Vi}) 
with a[ =Piai. 

Lemma 6.11 : Let (Gi,ai), (^2,02) and {G,a) satisfy eq.(58). Then 
a[G, a) = — (t(Gi, ai) + — (t(G2, 02)- 

"(G,a) a(G,a) 

Proof : Similar to that of Lemma 2.10. ■ 
In general, if (G, a) satisfies eq.(59) for some set of graphs {(Gj, Oj)}, we have, 

(^{G, a) = Y] d^Gi,ai)(^{Gi, tti). (60) 

"(G,a) ^ 

Theorem 6.12 : Every graph (G, a) having density matrix a{G, a) can be decomposed 

as in eq.(59), where a{Gi, ai) is a pure state. 

Proof : Same as that of Theorm 2.12. ■ 
Corollary 6.16 : A state of a /c-partite system is separable if and only if the graph 

{G, a) for a has the form 

(G, a) = [UiAf Hj^, (Gi, a>)] U [U,£ Hj^, (Gj, aj)] 

Ur^£[U,£Hj^,(Gj,ai)] (61) 

Proof : Same as of Corollary 4.11, where we refer to Theorem 6.4 instead of Theorem 
4.5 and Lemma 6.11 instead of Lemma 2.10 and eq.(60) instead of Remark 2.11. ■ 

Corolary 6.16 says that Werner's definition [1] of a separable state in C^^ ^C^^ (8)C^^ 
•••(g) C^*^ system, can be expressed using corresponding graphs. 



6.4 Representation of a hermitian operator (observable) by a 
graph 

In order to represent a general hermitian matrix corresponding to a quantum observable 
A we lift the requirement that the Laplacian be positive semidefinite and Tr[74] = 1. In 
other words we take the generalized Laplacian as the matrix for the graph. 

Given a Hermitian matrix A, the algorithm 6.1 can be implemented to get its graph 
(G, a). The corresponding observable A of a graph (G, a) is 

m ^ s 

k=i v2 
Example (12) : Give the graph of ax and ay. 
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(1) 



1 

1 



The corresponding graph of is shown in Figure 19a 

, . _ r -i 1 _ r e-'^/2 
^'""y ~\% J ~ [ e'^l^ 
The corresponding graph of Oy is shown in Figure 19b 



(a) 

Figure 19 



(b) 



Using Equation (46) to get the operators from graphs 



a, = 2P[^(|1) - e-/^|2))] - P[|l)] - P[|2)] = -z|l)(2| + ^|2)(1| 



1 

1 



-i 

1 



7 Some graphical criteria for the positive semidefi- 
niteness of the associated Laplacian 

In this section wc address the following question. Given a graph, can the positive scmidcfi- 
niteness of the associated Laplacian be determined using the topological properties of the 
graph? A general answer to this question seems to be difficult because the theory of 
weighted graphs, with negative and complex weights is almost unavailable. Many results 
obtained for simple graphs do not apply to the weighted graphs with real or complex 
weights. Nevertheless, we give here the above mentioned criteria in some special cases. 

Lemma 7.1 : Let {G,a) be a graph with real or complex weights, having one or 
more nonisolated vertices with degree zero. Then the Laplacian of {G, a) is not positive 
semidefinite. 

Proof : Such a graph (G, a) has one or more zeroes along the diagonal of its Laplacian 
with nonzero entries in the corresponding rows. However, a hermitian matrix with one or 
more zeros in its diagonal has at least one negative eigenvalue unless all the elements in 
the corresponding rows and columns are zero [18]. ■ 
Lemma 7.2 : Let (G, a) be a n vertex graph with real weights, having at least one 
loop and let the weights on all the loops be negative. Then the Laplacian of (G, a) is not 
positive semidefinite. 
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Proof : For the given [G, a) and some x in we have 

k t 

where the first sum is over edges and the second sum is over loops. It is easy to check 
that x'^[Q{G, a)]x < for x = (1 1 1 ... 1)^. ■ 
Lemma 7.3 : Let (G, a) be a graph without loops satisfying a{u, v) — au^e"^"", ^ 
27m). Then the associated Laplacian is positive semidefinite. 

Proof : This follows directly from Theorem 6.4. ■ 
Observation 7.4 : Let (G, a) be a graph satisfying a{u, v) = a,j„e*'^"" and a{{v, v}) > 
for all vertices in V{G,a). Then the associated Laplacian is positive semidefinite. 
Proof : The Laplacian is a hermitian matrix which is diagonally dominant. Therefore, 
by Gersgoin circle criterian [14, 15, 19] it is positive semidefinite ■ 
On a n vertex graph (G, a), we define a new graph operator Q{ui) which deletes the 
vertex Ui and rolls the edges incident on Ui into loops with same weights on the edges 
connecting neighbors of Ui as shown in Figure 20. We call the resulting subgraph principal 
subgraph. The Laplacian of the principal subgraph obtained by operating ©(itj) on {G, a) 
is the principal submatrix of the Laplacian of {G, a) obtained by deleting ith row and ith 
column. 




Figure 20 



Lemma 7.5 : If one or more principal subgraphs of (G, a) are not positive semidefi- 
nite, then (G, a) is not positive semidefinite. 

Proof : This follows from the result that all the principal submatrices of a positive 
semidefinite matrix are positive semidefinite[15]. ■ 
Lemma 7.6 : Let (G, a) be either a n vertex tree (n > 2) or a n vertex cycle (n > 4). 
We assume that there are no loops in (G, a) and that a{u, v) is real for all {u, v} e E{G, a) . 
Then (G, a) has a positive semidefinite Laplacian if and only if a{{u,v}) > for all 
{u,v) e E{G,a). 

Proof : Only if part : We prove that a{{u,v}) < for any one {u,v} G E(G.a) =^ 
Q{G, a) ^ 0. Let (T, a) be a tree with vi, . . . ,Vn vertices, and let {vi, Vi+i} be an edge in 
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(T, a) with negative weight a({fj,f j+i}) < 0. We operate on (T, a) by 0(fj+i). There are 
two possibiUties. If f j+i is a leaf , we get only one component with a negative weighted loop 
on Vi- By Lemma 7.2, the Laplacian of this principal subgraph is not positive semidefinite 
and by Lemma 7.5 the Laplacian of (T, a) is also not positive semidefinite. If Vi+i is not a 
leaf then 6(uj+i) will result in two or more principal subgraphs. The principal subgraph 
containing vertex vi is a graph having one loop with negative weight. By Lemma 7.2 
the Laplacian of this principal subgraph is not positive semidefinite and from Lemma 7.5 
Q{T,a) ^ 0. Let Cn be an n-cycle and let a{{vi,Vi+i\) = a < 0. We operate by Q{vi+i) 
which results in a n — 1 vertex path , say Pn-i with vi having a negative loop and 
having a positive or negative loop. If both the loops are negative we can use Lemma 7.2 
and Lemma 7.5 in succession to show that (5(C„,a) ^ 0. Suppose the loop on t'j_|_2 is 
positive . Then for some x e R^~^ we have 

n-2 

x^Q{Pn-i, a)x = ^ ai^j^{xi^ - Xj^f + a{vi+i, Vi+2)xl.^^ - \a\xl. 
k=i 

It is straightforward to check that x^Q{Pn_i, a)x < for = (1 1 ...0 1 ... 1), 
that is a vector x with all components 1 except Vi^2th component which is zero. Thus 
Q{Pn-i,a') ^ 0. By Lemma 7.5 Q{G,a) ^ 0. 

If p£irt : Assume Q{G, a) ^ 0. This imphes that there exists at least one x E 
satisfying 

x^(5(G,a)x = ^a(ife, - Xj^^f < 

k 

Since (xj^ ~Xjk)'^ ^ for the above inequality is satisfied only when a{ik,jk) < 
for some k. This proves the if part. ■ 

We observe that the proof of if part applies to all graphs as it should. 

Lemma 7.7 : Let all loops on a graph (G, a) have real positive weights. Let every 
edge {u,v} e E{G,a) having a{u,v) < be common to pair of C3. Let all such pairs of 
C3, each containing a negative edge be disjoint. Let all the edges in each pair of C3, other 
than the contained negative edge have positive weights satisfying a{u, v) greater than the 
absolute value of the weight on the negative edge. Then the Laplacian of (G, a) is positive 
semidefinite. 

Proof : Consider a negative edge common to two Cs's as shown in the Figure 21. 

By hypothesis bi > a,i = 1, 2, 3, 4. We can write 6j = a + q, q > 0, i = 1, 2, 3, 4. We 
can decompose this graph as the edge union as shown in Figure 22. 

The first graph on RHS has all positive weights and hence has a positive semidefinite 
Laplacian. It is straighforword to check, that the second and the third graphs on RHS 
correspond to the projectors P[-i=(|j) — 2\l) + \k))] and P[-l=(|j) — 2\i) + \k))] respectively. 
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Hence they have positive semidefinite Laplacians. The Lapalcian of the graph on LHS 
is the the sum of the Laplacian of the graphs on RHS (Lemma 2.10), each of which is 
positive semidefinite. But we Know that the sum of positive semidefinite matrices is a 
positive semidefinite matrix [15]. Now the graph {G,a) can be written as edge union of 
the factors (spanning subgraph) as Figure 21 (possibly more than once) and the remaining 
factor which has all positive weights. The Laplacian of each factor is positive semidefinite 
and the Laplacian of the given graph, being the sum of positive semidefinite matrices, is 
positive semidefinite. ■ 
Lemma 7.8 : If all the negative edges of a real weighted graph {G, a) occur as in the 
following subgraph as shown in Figure (23), where Cj > 6 ; i = 1, 2, . . . , 8 and 6 > a > 
then the associated Laplacian is positive semidefinite. 

Proof : We can decompose the above graph into factors as shown in Figure(24) 

From the graphical equation in Figure (24) we see that the first factor on RHS corre- 
sponds to -P[|— )|— )], second factor has a positive semidefinite Laplacian from Lemma 7.7 
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b > a > 

Figure 23 




and the third factor has a positive semidcfinitc Laplacian as it has all positive weights. 
Since this graph occurs (once or more) as disjoint subgraphs of {G, a) it can be written 
as edge union of one or more of these subgraphs and the remaining graph containing only 
positive or complex edges. Since each of these has a positive semidefinite Laplacian, [G, a) 
also has a positive semidefinite Laplacian. ■ 
Lemma 7.9 : Let {G'^",a) be a complete signed graph with weight function aij e 
{—1,1} without loops on 2" vertices n > 1. Let Ei denote the set of edges incident on 
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ith vertex {\Ei\ = 2" — 1) and let E^,Ej^ denote the sets of edges incident on the the 
ith vertex with weight +1 and —1 respectively, {Ei = Ef + E^). Let (G^^.a) satisfy 
the following condition (i) \E^\ = 2"~^ — 1, i = 1, 2, . . . ,2", so that the degree of every 
vertex= 1. Then (G^", a) corresponds to a pure state in 2" dimensional Hilbert space. 
Proof : We need to prove that condition (i) in the statement of the Lemma can be reahzed 
for all n and that the resulting signed graph corresponds to a pure state for all n. We 
use induction on n. It is clear that the assertion is true for n = 1 with the corresponding 
pure state given by P[;^(|l) — |2))]. Now assume that assertion (that is condition (i) 
and purity of the corresponding state) is true for n = k. For the graph corresponding to 
n — k + 1 with a) \ — 2^+^ consider the modified tensor product 

{G^,a)B{G^\a) = {G^'^\a) ^ {C{G^,a) ^ /:r]{G^\a)} + {jC{G^,a) ^ M{G^\a)} 

+{A^(G'^ a) ® C{G^\a)} + {^(G'^ a) ® n{G^\ a)} (59) 

where C, A/" and Q arc graph operators defined in equation (20b) and G^, G'^^ arc graphs 
with number of vertices 2 and 2^ respectively. Note that the last term corresponds to 
an empty graph as G does not have any loops. Since the modified tensor product of 
two complete graphs is also a complete graph, (G'2'=+^ , a) is a complete graph. Therefore 
\Ei{G'^''^\a)\ = 2'=+! - l,i = 1, 2, . . . , 2^=+^ To show that condition (i) is realized for 
^^2*^+1 ^ given the induction hypothesis, we note that the first term in equation (59) 
contributes \Ef{G'^'',a)\ negative edges to (l,i)th vertex in [G'^'^^^) and the third term 
contributes |£'j~(G^'°, a)| negative edges, while the other two terms have no contribution. 
Therefore 

|Er.(G^'"\a)| = \Et{G'\a)\ + \Er{G^\a)\ =2^-1 

Similarly, the first three terms contribute |£Jj~(G^'°, a)|, 1 and \E^{G'^^ ^a)\ positive edges 
to (1, i)th vertex. Therefore, 

\Et^G-''^\a)\ = \E-{G''\a)\ + \E+{G''\a)\ + 1 = 2^'^ - 1 + 2'^-' + 1 = 2^= 

and similarly for i?2i- That {G^'"^^ , a) corresponds to pure state follows from the fact that 
the state corresponding to the modified tensor product of two graphs is the tensor product 
of the states corresponding to the factors. Since the state for {G^'° , a) is pure by induction 
hypothesis and {G'^,a) is pure, the preceding statement means that (G'^'°^\ a) is a pure 
product state. ■ 



8 Summary and Comments 

Following is a brief summary of the main features of the paper 
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(i) We have given rules to associate a graph to a quantum state and a quantum 
state to a graph, with positive semidefinite generahzed Laplacian, for states 
in real as well as complex Hilbert space(2.1 and 6.1). 

(ii) We have shown that projectors involving states in the standard basis are 
associated with the edges of the graph (Theorems 2.7 and 6.4) 

(iii) We have given graphical criteria for a state being pure. In particular we 
have shown that a pure state must have a graph which is a clique plus 
isolated vertices (Theorem 2.3, 2.4, Remark 6.2) 

(iv) For states in a real Hilbert space , we have given an algorithm to construct 
graph corresponding to a convex combination of density matrices, in terms 
of the graphs of these matrices (2.3.2). 

(v) We have defined a modified tensor product of two graphs in terms of the 
graph operators £, 77, A/", Q and obtained the properties of these operators 
(4.2, 6.2). We have shown that this product is associative and distributive 
with respect to the disjoint edge union of graphs (corollary 4.7, Remark 
6.11). 

(vi) We have proved that the density matrix of the modified tensor product 
of two graphs is the tensor product matrices of the factors. (Theorem 4.5, 
6.10 ). For real density matrices, we show that a convex combination of 
the products of density matrices has a graph which is the edge union of the 
modified tensor products of the graphs for these matrices (corollary 4.11). 
Thus we can code werner's definition of separability in terms of graphs. 

(vii) We have generalized the separability criterion given by S. L. Braunstein, 
S. Ghosh, T. Mansour, S. Severini, R.C. Wilson [2] to the real density ma- 
trices having graphs without loops Lemma (4.16). 

(viii) We have found the quantum superoperators corresponding to the basic 
operations on graphs, namely addition and deletion of edges and vertices, 
it is straightforward to see that all quantum operations on states result 
in the addition / deletion of edges and/ or vertices , or redistribution of 
weights. However, addition / deletion of edges / vertices correspond to 
quantum operations which are irreversible, in general. Hence it seems to 
be difficult to encode a unitary operator, which has to be reversible, in 
terms of the operations on graphs. Further , graphs do not offer much 
advantage for quantum operations which only redistribute the weights, 
without changing the topology of the graph, as in this case the graph is 
nothing more than a clumsy way of writing the density matrix. 
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(ix) In section 6, we generalize the results obtained in sections 2 - 4, to quantum 
states in a complex Hilbert space, that is, to the density matrices with 
complex off-diagonal elements. In fact, all the results previous to section 
5 go over to the complex case, except Lemma 4.16. Many of these results 
have been explicitly dealt with (e.g. Theorem 6.4, Remark 6.2, section 6.2 
etc) . We have also given rules to associate a graph to a general hermitian 
operator. We believe that any further advance in the theory reported in 
this paper will prominently involve graph operators and graphs associated 
with operators. 

(x) Finally, we have given several graphical criteria for the positive semidefi- 

niteness of the generalized Laplacian associated with a graph. Note that 
by Lemma 7.3 and observation 7.4 all graphs with complex weights, either 
without loops or with positive weighted loops have positive semidefinite 
generalized Laplacians. This characterizes a large class of graphs coding 
quantum states. 

This paper is essentially a generalization of the work by Braunstein, Ghosh and Sev- 
erini [7] in which the idea of coding quantum mechanics of multipartite quantum systems 
in terms of graphs was implemented. The motivation in both Braunstein, Ghosh, Severini 
and this paper is to explore the possibility of facilitating the understanding of mulipartite 
and mixed state bipartite entanglement using graphs and various operations on them. 
Whether such a goal can be reached is too early to say. In order to code arbitrary quan- 
tum states and observables in terms of graphs, we have to deal with weighted graphs with 
real or complex weights. Unfortunately, a mathematical theory of such graphs is lack- 
ing. Many results pertaining to simple graphs are not available for such weighted graphs. 
We hope that, through the need of understanding entanglement and related issues the 
mathematical structure of weighted graphs gets richer and in turn gives a feedback to our 
understanding of entanglement. 
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